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Abstract 

The authors modify a non-perturbative approach based on the Principle of Minimal Sensitivity to 
calculate the periods of relativistic oscillators with even polynomial potentials. The optimization 
of the variational parameter is adapted by using the ultrarelativistic limit of the period as a 
boundary condition, which determines the values of additional free parameters introduced. General 
approximations for the potentials \ + Yl™=i is^T an d prove to be accurate over the whole 
solution domain and even for large values of m. 
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I. INTRODUCTION 



Interest in the relativistic 



larmonic oscillator has either focused on its mathematical and 



physical characteristics 



2l, 14, [5|, la 13] or on analytic approaches of solving its non-linear 



equation of motion, given non-dimensionally as 



X + x(l — X )a = 0. 



(1) 



The latter is motivated by the general significance of the study of strongly non-linear oscil 
lators, particularly in developing efficient approximation methods to treat them 



Recently, the non-perturbative methods of harmonic balance jll, 12, is| and homotopy 
[bl [loT ] have been succesfully applied to give analytic approximations of the oscillation period 
and periodic solution of EqJTJ Unlike traditional perturbation theory, these asymptotic 
techniques do not rely on any linear term or physically small parameter. Hence their results 
are unrestricted by a radius of convergence. For both methods, the approximations had 
minimal error over the entire domain of the particle's maximum velocity, x max G [0, 1) 22]. 
However, the applicability of these approaches has thus far been restricted to the linear 
potential. 

In this work we contribute a novel method for computing analytic approximations of the 
period of strongly anharmonic oscillations described by 



x + —(l 
dx 



0, 



(2) 



for an even polynomial potential (p(x) of general form. We adopt a non-perturbative ap- 
proach that has been developed and applied by Amore et al. in computing periods of 
non-linear oscillators 



16, 



171 ] and formally equivalent systems such as the deflection angle 



in gravitational lensing [18]. The method operates by converting the relevant integral into a 
series carrying a variational parameter u, which can be chosen using the Principle of Mini- 
mal Sensitivity (PMS) [li] ]. In our case the standard procedure for choosing u) is modified 
by treating the existing harmonics as free parameters. These are then assigned values by 
imposing the ultrarelativistic limit of the period as a boundary condition. Although the 
necessary computations are simple, the resulting approximations are compact and accurate 
over the whole solution domain. 
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II. METHOD 



As applied to the relativistic harmonic oscillator the qualitative analysis of ordinary 



differential equations detailed in Ref. 



20j can be used to show that the solution to Eqf2] is 



indeed periodic. That is, the phase space trajectories given by 

dy dd> (1 — y 2 )^ 

-T + -T- — = 0, 3 

dx dx y 

where y = x, are closed and have the strip structure: 

— oo < x < +oo and - 1 < y < +1. (4) 

The period of this system can be derived from the generalized energy-momentum relation 

E = y/p* + l + <j>(x), ( 5 ) 

where E is the total energy and p is the relativistic momentum. In all our equations the 
rest mass and speed of light are set to unity. 

By noting that the total energy is the sum of the potential energy at the amplitude A 
and the rest mass, i.e. E = 4>(A) + 1, it can be shown from Eqj5]that 

<P(A) - cj>{x) + 1 = (1 - y 2 yK (6) 

which leads to the following expression for the period: 

T =2 f A dx * (A) ~ M + 1 (7) 

S-A ^{4>(A)-^(x)]{4,(A)-^( x )+2\ 

Since the period is a function of the oscillation amplitude, it would be useful to express 
the physical regimes described by x max in terms of this parameter. Setting x = in Eqj6] 
gives 

l + HA) ' [) 

This relation maps the domain x max £ [0, 1) unto A G [0, oo). 

For an even potential (fi(x), the period from EqJ7]can be written as 



1 dx 



T = 2 =. (9) 

-A ^R(x)[A 2 -x 2 ] 
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We approximate R(x)(A 2 



x- 



with the solvable harmonic potential oj 2 (A 



x 



manner that is essentially a Linear Delta Expansion 

rl dx 
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2l|: 



T 



m a 



(10) 



-A ^[A 2 - x 2 }[u 2 + 5(R(x) -u 2 )] 

In the above equation, when 5 = we have the case of simple harmonic motion, while for 
5 = 1 the system reduces to the original potential of EqJHl 

With a change of variable x = AcosO, the above equation becomes 



T 



2 



de 



(ii) 



where the substitution A(9, ui) = —1 + has been made. A binomial expansion of the 
integrand in powers of A(9,u) then gives 



r =^EU) r* n (0,")d9, (12) 

where we have set 5=1. 

Evaluating the above series to a finite order will leave an explicit w-dependence that 
should be eliminated. Indeed this is the basis of the PMS, requiring that |^ = 0. As shown 



o 



in Ref . [17| , the condition is equivalent to choosing u such that 

A n (6,uj)d6 = (13) 
for an nth order optimization. To first order, the above criterion is equivalently 

VVW), (14) 



where (R(0)) denotes the mean value of R(9). 

For classical oscillators Eq{14]is implemented easily, yielding remarkable results [17|. In 
the case of relativistic oscillators however, the expressions not only become cumbersome 
but are divergent at large oscillation amplitudes even with higher order computations. For 
example, we find that the resulting uj computed for the relativistic quartic oscillator ap- 
proaches unity as A —>■ oo. This marks an ineffective variation since the choice to = 1 gives 
the exact untreated series. Hence for any order of calculation, the error of approximation 
increases rapidly beyond Ami [x max ~ 0.82). 



The alternative criterion given by Eq.36 of Ref. [17|] is also found to be problematic because 
the resulting u goes to infinity as A — > oo. Since all terms in Eq{12]are proportional to powers 
of -, then any finite order approximation incorrectly goes to zero in the ultrarelativistic limit 
A -> oo. 

The correct large amplitude behavior can be verified by noting from Eqj6] that 4>(A) — 
4>(x) >> 1 in this regime. Hence, EqJT] simplifies to what is expected for a photon bouncing 
in a box of width 2A: 

T ~ 4A (15) 

The preceding limit will serve as the asymptotic boundary condition that determines our 
variational parameter. 

We begin constructing our solution by noting that the zeroth order truncation of the 
series in EqJT2l is 

T = X (16) 

to 

We then assume the form of uj that directly satisfies Eq JT3l namely 



y/Wj. (17) 



In the above equation uj is still dependent on the variable of integration 9 since harmonics 
of the form cos 2fc 9 are present. To remove this, we treat the harmonics as free parameters, 
denoting cos 2k 9 as X 2 k- 

Next, we look at the limit of u as A — > oo. It can be shown that for an even potential 
with leading order x 2m , the large amplitude behavior is 



„ _ . 1 Y,n=o(^2n ~ A 2(m+n) ) ^ 
A y 1 — 2\2 m + ^4m 

Substituting EqsJT5] and [18] in the limit A — » oo of Eq|16] will yield a relation between all 
orders of free parameters that exist in R. This relation can be used recursively to provide m 
equations. However for a non-linear potential, the system of equations is under-determined 
since there are a total of 2m free parameters. We remedy this problem by factoring out 
1 — ^2m from the numerator and denominator in Eqfl8| giving us 



w ~iVt^t- (19) 
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Thus the large amplitude limit of Eq{16] yields the relation 

A 2m = 1 - — 2 A 2 n- (20) 

n=0 

Under the ansatz of setting the lowest order free parameter to be its trigonometric average, 
i.e. A 2 = \, the m = 2 case of the above equation determines A4. These values are then 
used in the m = 3 equation to solve for A6 and so on, effectively building the whole set of 
m free parameters from the bottom up. This iterative computation yields that for integers 
n > 2 the free parameters are given by 

7T — 6 / 7T — 4 N 



^ = —4{— ) ■ < 21 > 

We use these values in R and determine u by Eq{T71 Then by Eq{16]the resulting approxi- 
mation is 

T= (22) 

V R(^2, A4, A2 m ) 



III. APPLICATIONS 



As a first example we apply the method to the linear potential. Plugging in <fi(x) 
in EqJT] and following the general derivation previously outlined yields 

R = ± + A\l-\ 2 ) 

4 + 4A 2 (l-A 2 ) + A 4 (l-2A2 + A 4 )- 1 ; 

In calculating the free parameters present in R we utilize the relation given by Eq{18] since 
for m — 1 the system is not under-determined. Accordingly, the ultrarelativistic limit of 
Eq HE] gives the relation 



WW- 

By using A2 = \ in the above equation, we determine A4 = We then substitute these 
values into Eqj23] and use Eqj22l giving us the compact expression for the period 

A^/A± + (2 + gj g 

Th ~ TfT^ — • (25) 
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FIG. 1: (Color online) Relative error of Eqj25] with respect to the exact result given by Eql261 
Shown for comparison are the relative error of the second and third order computations of the 
harmonic balance method [12I . and the first to third order computations of the homotopy 
method [15I ]. 



We compare this with the exact period given by 



4 (4 + A 2 ) E 



4+A 2 



- 8K 



i+A 2 



exact 



(26) 



where E[-] and K[-] are the Elliptic Integrals of the First and Second kind, respectively. 



Figure [T 
balance 



shows the relative error of Eqj25] along with those of the results of harmonic 



131 ] and homotopy [15j]. Despite the simplicity of our calculation, the third 



order approximations of both these methods are only more accurate by 0.1 — 0.2% over a 
finite range of A. Our solution is the only approximation whose error goes to zero at the 
ultrarelativistic limit. An advantage is also apparent in the resulting expression itself, that it 
is compact and mathematically elementary, while those of harmonic balance and homotopy 
are lengthy and even require the evaluation of Elliptic Integrals. 

The most note-worthy point of our method is its extensive applicability. For our second 
example we calculate a general approximation for the periods of relativistic oscillators with 
the potential (j)(x) = + t^-. As we have done previously, we obtain the function R by 



following the derivation starting from EqJTl The result is 



m— 1 



n=0 V 17 

where 

Zi = mA 2 (l - A 2 ) + A 2m (l - A 2m ). (28) 

Note that this expression has been factored in a manner compatible with Eq{19] so that the 
existing free parameters are defined by Eqf2~ll Then the approximation simply follows from 
Eqj22l Explicitly, the result is 

T a = 2n(2m +Zl ) (jjg) 

yj (4m + Z x ){m + A^\\ + ^(1 - (4^) m ~ 2 )]} ' 

where we have evaluated the geometric sum appearing in Eqj27l To give specific examples 
of this general result, the approximations for the periods of the relativistic quartic [m = 2) 
and sextic (m = 3) oscillators are respectively 



2V2[6A 4 + (4 + A 



2 % 2 



TquarUc ^(4 + 3A 2 )[QA* + (8 + A 2 )7r 2 ] (30) 



and 



24?r 4 + GA 2 n 4 + A 6 (40tt 2 - 96) 

sextic ~ ^[6tt 2 + A\hn 2 - 12)][24tt 4 + 3A% 4 + A 6 (20vr 2 - 48)] ' ' ' * ' ' 

The approximations for the cases m = 2,3,4,20,500 and the harmonic potential are 
fitted with the numerically integrated exact results in Figure [2 The relative errors shown 
in the inset are remarkably low even for the highly non-linear case of m — 500. For such 
large values of m the period exhibits a sharp jump from what is characteristically harmonic 
at amplitudes A < 1 to the ultrarelativistic limit for A > 1. 

As a third example, we derive a general approximation for the periods of relativistic 
oscillators with the potential (j)(x) = Y^n=i 1r~ - Again, the crux of the calculation is deriving 
a general form of R and factoring it compatible with Eq{19] so that the free parameters are 
determined by EqEU Then by using Eql22l the following result is obtained 

2 ^ (1 + Z2) (32) 

;2 + z 2 ){\ + e: =2 ^[§ + ^(i - (^r- 2 )}} 
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FIG. 2: (Color online) Numerically integrated exact results (approximations) for the cases m = 
2, 3, 4, 20, 500 and the harmonic potential are denoted by E m (T m ) and Eh (T#), respectively. The 
inset shows relative errors labeled by the respective subscripts. 

where 

^E-^ 1 -^)- (33) 

n=l 

For m = 2 this reduces to Eql30l the period of a relativistic quartic oscillator. For m — 3, 
it provides the following approximation to the period of a relativistic particle oscillating 
under the potential cj)(x) = ^ + x + T 

_ 2 v / 2[-48A 6 + 2A 4 (9 + 10yi 2 )7r 2 + 3(4 + yl 2 )7r 4 ] 

2,4,6 ~ v^-24,4 4 + (12 + 9A 2 + 10A 4 )^ 2 ][-48A 6 + 2A 4 (9 + l(L4 2 )7r 2 + 3(8 +^4 5 )i i J" 

(34) 

Figure [3] shows the approximations for the cases m = 3,5,10 and 100, fitted with the 
numerically integrated exact results. Notably, the inset shows low relative errors whose 
maximum do not increase with increasing m. 

Lastly, we treat relativistic oscillators with the intrinsically non-linear potential y~- 
Using the same procedure of deriving a general expression for R and using Eqsj21] and [22l 
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FIG. 3: (Color online) Numerically integrated exact results and approximations for the cases 
m = 3,5, 10 and 100 are denoted by E m and T m , respectively. The inset shows relative errors 
labeled by the respective subscripts. 



Tc = 27r[2m + A 2m (l-A 2m )] 



the result is 



A— y [4m + A^(l - A 2m )]{§ + ^[1 - (^) m - 2 ]} 
For example, setting m = 2 in the equation above gives the following approximation for the 

4 

period of a relativistic oscillator with potential \ 

= 12A^ + ^ 
V9AH 2 + 12A% 4 

The approximations for the cases m = 2,3,4 and 20 and their numerically integrated 
exact results are shown in Figure HI The method is able to capture the infinities at both 
regimes of the oscillation period. 

IV. CONCLUSIONS 



Our modification of a PMS-based non-perturbative method offers an economical proce- 
dure of calculating compact and accurate analytic expressions for the period of an oscillating 
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FIG. 4: (Color online) Numerically integrated exact results and approximations for the cases 
m = 2, 3, 4 and 20 are denoted by E m and T m , respectively. The inset shows relative errors labeled 
by the respective subscripts. 

relativistic particle. It is able to attack even polynomial potentials of general form, including 
those that are intrinsically non-linear. The mechanics of the method itself guarantees the 
accuracy of the solution even at the ultrarelativistic limit. 

J. P. E. acknowledges support from a National Institute of Physics Faculty Grant. 
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